A surface electrode point Paul trap 
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We present a model as well as experimental results for a surface electrode radio-frequency Paul 
trap that has a circular electrode geometry well-suited for trapping of single ions and two-dimensional 
planar ion crystals. The trap design is compatible with microfabrication and offers a simple method 
by which the height of the trapped ions above the surface may be changed in situ. We demonstrate 
trapping of single 88 Sr + ions over an ion height range of 200-1000^m for several hours under Doppler 
laser cooling, and use these to characterize the trap, finding good agreement with our model. 



I. INTRODUCTION 

Radiofrequency (rf) traps have been applied exten- 
sively in a large variety of scientific studies over the past 
six decades. Originating from mass spectrometry [![, 
they have then been applied in fields such as metrol- 
ogy 0! quantum information science [!, H[ and cold 
molecular physics [5, 6], to mention but a few. 

Traditionally, such devices have been rather bulky, 
three-dimensional structures that required precise ma- 
chining and careful assembly. Recently, however, the 
four-rod linear Paul trap (see Fig. [1^) has been trans- 
formed into a two-dimensional structure, with all elec- 
trodes in a single plane above which ions can be 
trapped @- This new class of so-called surface traps 
offer a tremendous advantage over their predecessors in 
that electrodes can be defined lithographically with ex- 
tremely high precision and that construction can leverage 
the techniques of microfabrication, with the possibility 
of incorporating the technology of CMOS for integrated 
control hardware [1, H[ . These aspects are particularly 
attractive to applications in quantum information pro- 
cessing where limitations are currently, by a large degree, 
pertaining to the scalability of devices for trapping as well 
as certain elements of infrastructure such as optics, laser 




FIG. 1. Comparison of the traditional four-rod linear Paul 
trap (a) to the point Paul trap (b). The latter achieves 
quadrupole ion confinement through rf on a single, ring- 
shaped electrode. Dashed lines suggest how cylindrical ele- 
ments - such as optical fibers - may be introduced to the 
point Paul geometry. 
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light delivery, and control electronics. 

In this paper we study a novel type of rf surface trap 
with a high degree of symmetry in its electrode geometry 
The generic geometry of this trap, which we shall refer 
to as the point Paul trap, is shown in Fig. QJ>, and may 
consist of any number of concentric electrodes of arbi- 
trary widths to which different voltages can be applied. 
This design originated in a study of surface electrode 
traps fiol. TlTj] . but was subsequently strongly inspired by 
work on planar Penning traps [l2| , where a similar geom- 
etry was used to create a static electric quadrupole field 
that, when combined with a strong homogeneous mag- 
netic field, gave rise to a confining potential above the 
surface of the electrodes. The point Paul trap also bears 
close resemblance to the rf ring- and the rf hole-traps [13| ; 
however differs in that the ion is trapped above the sur- 
face of the electrodes as opposed to in-between, which 
makes this geometry better suited for microfabrication. 

A consequence of the azimuthal symmetry of the 
electrodes is that the rf-field exhibits a nodal point 
rather that a nodal line as in the linear Paul trap 
and that the confining fields originate exclusively from 
the rf-potential, rendering the addition of dc-potentials 
nonessential for anything but the compensation of stray 
charges on the trap. This makes the point Paul trap 
well-suited for confinement of single ions, which may then 
reside at the rf-nodal point where the amplitude of the 
rapidly oscillating rf-ficld vanishes. 

The ability to fabricate these traps in a scalable fash- 
ion makes them attractive for realizing large arrays of 
single ions in independent traps that may be utilized for 
a quantum processor, provided the individual ions can be 
interconnected, e.g., through optical fibers [3 EH- On 
this aspect, the axial symmetry of the trap lends itself 
well to integration of such fibers and potentially other 
optical elements that also possess axial symmetry. The 
fiber, for instance, may be introduced through the elec- 
trodes directly beneath the ions with minimal perturba- 
tion of the trapping fields. 

Another possible application of this trap is in the 
field of quantum simulation. While classical computers 
are unable to efficiently simulate coupled spin systems, 
such simulations may be implemented using a quan- 
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turn mechanical system of effective spins, such as a two- 
dimensional lattice of interacting ions. The resulting po- 
tential of the point Paul trap provides ion crystals with 
exactly the requisite two-dimensional planar structure. 
As such, the system could be used to simulate e.g. a 
frustrated spin system [l6|, [ItJ , as was demonstrated re- 
cently [HI]. 

We also find that our trap design is ideally suited for 
realizing a scheme by which the height of a single trapped 
ion above the trap surface is varied in-situ. This capa- 
bility may prove extremely useful in the search for the 
origin of anomalous heating in ion traps - a problem 
currently impeding the advancement of quantum com- 
putation with trapped ions 0, Ho| . It also provides a 
general technique by which oven contamination of the 
trap can be minimized by loading further away from the 
trap surface, and subsequently bringing the ion to the 
desired trap height. 

This paper is organized as follows: in Section |TT] we 
present a model for the planar point Paul trap, derive 
analytic expressions for the relevant trapping parame- 
ters, present full numerical results of trap optimization, 
and consider a scheme for the variation of the ion height 
above the trap surface. In Section UTT1 we describe our ex- 
perimental setup used for the verification of the model; 
and in Section IIVI we present experimental results for 
trapping of single and few ions in a printed circuit board 
(PCB) implementation of the point Paul trap. 

II. POINT PAUL TRAP MODEL 

We proceed with a general treatment applicable to an 
arbitrary number of circular electrodes, and then focus on 
a particular geometry that we will study experimentally 
later in Sections IIIII and IIVI At the end of this section 
we also describe a scheme for variation of the ion height 
above the trap surface. 

A. Potential from annular planar electrodes 

We begin with the general solution to the Laplace 
equation in charge-free space and ex pre ss this in cylin- 
drical coordinates for z > 0, yielding [21| 

°° />oo 

$( 2 ,p,^) = V / J m (kp)e- kz (1) 

m=0 J ° 

x [A m (k) cos (m<j)) + B m (k) sin (m<p)] dk. 

where J m (kp) are the usual Bessel functions and A m (k) 
and B m {k) are coefficients to be determined based on 
the boundary conditions of the problem. Based on the 
azimuthal symmetry of Fig. [2J Eq. (TTJ) further simplifies 

to du 

p OO 

Hz,p)= .h{kp)e- kz A {k)dk. (2) 
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FIG. 2. The generic layout of the point Paul trap, which con- 
sists of concentric annular electrodes with arbitrary widths, 
illustrated in a cylindrical coordinate system. 



In turn, Aq can be expressed as Ao(k) = J^" Ai(k), 
where 

Ai{k) = k / pJ (kp)Vi(p)dp. (3) 

J Oti 

All information about the electrode geometry is now in- 
cluded in the Ao coefficient, which we have in turn writ- 
ten as a sum of n sub-coefficients, each accounting for the 
effect of a single annular electrode i with inner radius cti, 
outer radius (3i and a voltage Vi , which we shall assume is 
constant across the electrode. The integral of Eq. ([3]) can 
be evaluated using the identity for the Bessel functions 
Jo vJo(v)dv = uJ\(u) to give 

Mk) = Vi [ftJi(feft) - otiJiikcti)] ■ (4) 

This completes the general treatment of the problem: 
The electric potential above a surface at z = with 
n concentric circular electrodes, each with independent 
voltages Vi and inner and outer radii of on and Pi respec- 
tively, is given by Eq. ([2]) with Ao(k) = Y^i M{k), where 
the Ai coefficients are given by Eq. (fj|. 



B. The three-electrode point Paul trap 

While static potentials alone may provide confinement 
in one or two dimensions, Earnshaw's theorem dictates 
that this is accompanied by a defocusing effect in the 
orthogonal dimensions. In the work of Ref. [l2j], three- 
dimensional confinement was achieved via the addition of 
a static magnetic field to realize a planar Penning trap. 
Here, we use a time- varying rf field to achieve charge 
confinement as a planar Paul trap. Namely, we consider 
the simple geometry of only three electrodes defined by 
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the following boundary conditions: 

'0 for < p < a, 

$(z = 0, p) = { Vri cos (flrft) for a < p < b, 
for b < p < oo, 



(5) 



where V r { is the amplitude of the applied voltage and O r f 
is the frequency. The resulting potential then reads 



$(z, p, t) — V T { cos (fi, TI t)K(z, p), 



(6) 



where 



«(*,p) = / [6Ji(fc6) 
Jo 



aJ 1 (ka)}e^ kz J a (kp)dk. (7) 



In general, Eq. ([7]) has to be solved numerically. How- 
ever, for the case of p = the problem simplifies signif- 
icantly and an analytic solution can be obtained. From 
the symmetry of the problem we can infer that if a non- 
trivial field zero (E — 0) exists for z > 0, it will be located 
on the axis defined by p = and hence this scenario is 
worthy of attention. 

The on-axis potential is easily integrated to yield: 



k(z,0) = 



1 



1 



(8) 



Asserting that this potential provides for an electric field 
null at z — zq > 0, we expand <f>(z,t) around this point 
and use the resulting expansion to find the equation of 
motion for a particle of mass M and charge Q along the 
z-axis: 



m = 



q d$(z,t) 

M dz 
QV T 



(9) 
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- cos (firfi) x [f(a, b)(z - z ) + 0(z - z ) 2 



Provided \z — zq\ <C Zq, which is a reasonable assump- 
tion for a trapped ion, terms of second and higher order 
can be neglected, and the equation of motion takes the 
form of the well-known Mathieu equation. With proper 
rescaling of variables, Eq. © can be cast into the stan- 
dard Mathieu form |22| : 



z(t) + 2gcos(2r)5(r) = 0, 



(10) 



where we have substituted z = z — zq and r = fl rl t/2, 
and the Mathieu q-parameter has been defined as 



2QV lf 



(ii) 



Here, everything related to the trap geometry is collected 
into a single function /(a, b) of unit [length] -2 given by 



%{bi~aif{bi+aif 

/( a > b > = \ 71 a. , A 72—i 477- (I 2 ) 

y 03a3(t)3 +b3Q3 +a3j° 



Note that in this treatment a denotes the inner radius of 
the rf electrode and not the Mathieu a-parameter, com- 
monly used in the literature on Paul traps. The Math- 
ieu a-parameter, which corresponds to the inclusion of a 
dc-potential in the equation of motion (Eq. [9]), is ren- 
dered superfluous in the point Paul trap as full three- 
dimensional confinement is achieved in this geometry by 
the rf-field alone. 

When the trap is operated such that \q\ -C 1, the equa- 
tion of motion can be readily solved to yield 



Z(t) = (TO 



1 



| cos (firfi) 



cos (uj z t). 



(13) 



This is the usual result, familiar also from the four- rod 
linear Paul trap, where the motion is comprised of two 
distinct types of motion: a slow, so-called secular, motion 
with an amplitude ctq at a frequency 
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(14) 



and a superimposed, fast micromotion, with a lower am- 
plitude of ^qo-Q and at sideband frequencies of the rf drive 
ri r f . Neglecting the micromotion - a reasonable approx- 
imation for \q\ <C 1 - we can define an approximate har- 
monic potential that describes the ion motion near the 
quadrupole zero by the following 



■*{z)= l -Mul{z-z ) 2 



Q 2 Vr\ 

4MfE 



f(a,b)(z-z ) 2 , (15) 



thereby showing the charge-confining capabilities of the 
three-electrode point Paul trap, provided that f 2 > 0. 



C. Trap optimization and results 



While the harmonic potential of Eq. (|15[) provides an 
intuitive connection to the physical, time-averaged mo- 
tion of the trapped ion in the vicinity of the rf node, 
it does not reveal any information about the dynamics 
where the inequality \z — Zq\ <C zq is not satisfied. For 
instance, in a real device there is necessarily a finite trap- 
ping volume and trap depth. These quantities originate 
from the shape of the potential significantly beyond the 
harmonic region. In the limit q <C 1, the effective po- 
tential energy beyond the harmonic regime - commonly 
referred to as the pseudopotential - may be expressed di- 
rectly through the gradient of the electric potential, here 
written in terms of 



Q 2 Vr\ 

AMtl 2 



\Vk{z, P )Y 



(16) 



The solid line in Fig.[3]shows the pseudopotential for the 
planar point Paul trap. Superimposed is the harmonic 
approximation (dashed line) . Inserting the expression for 
k from Eq. ([7]) into the expression for the pseudopotential 
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FIG. 3. The harmonic potential approximation (Eq. 1 1 5 1 
dashed line) for the point Paul trap is shown on top of the 
pseudopotential (Eq. 1161 solid line). The latter identifies the 
trap turn-around point z max in addition to the trap location 
zo- Trap geometry (a,b) is chosen as in Table [I] 



(Eq. I16[) yields two physically meaningful extrema, 



zo = 



1 6 4 / 3 q 4 / 3 
6 2 / 3 + a 2 / 3 



(17) 



TABLE I. Results of numerical optimization of the trap depth 
for a fixed trap height zq. The choice of units allows for direct 
comparison with the three-dimensional linear Paul trap. 
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the four-rod linear Paul trap we have defined 
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2QV rf 



and £>4rod = 



rf^O 



(20) 



which corresponds to the g-parameter and the trap 
depth, respectively, for the three-dimensional four-rod 
Paul trap with an ion-electrode distance of Zq. 

From the optimization results of Table HI it is seen that 
the g-parameter of the point Paul trap is roughly a factor 
h of the four-rod linear Paul trap while the trap depth is 
about a factor 4^. By comparison, the surface electrode 
linear Paul trap design that has recently attracted much 
attention in the context of quantum computing (26l . 
has a g-parameter and a trap depth that are approxi- 
mately a factor -^j= and ^ , respectively, of the four- rod 

linear Paul trap |7j . 
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(18) 



Here, zq is the coordinate of the pseudopotential min- 
imum, and z m ax denotes the turning-point of the con- 
fining pseudopotential. The difference z max — zq can be 
taken as a linear measure of the effective trap volume. 
Furthermore, the corresponding trap depth can now be 
defined as D = *(z max ) - *(z ). Using Eq. ([TTJ) and 
Eq. (fT5)l , one finds that the trap depth is positive for all 
values of b > a > 0, and is equivalent to: 
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Q 2 v r 2 { 

4M0 2 f 
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b 2 (b 2 



z 2 )-' 

max/ 



(19) 

The trap depth is a reasonable quantity to be opti- 
mized in the design of the point Paul trap. However, 
unconstrained optimization of Eq. (|19p over (a, b) will 
influence not only the trap depth but also the trap height 
Zo above the surface through Eq. (fTTj) . Often in experi- 
ments, the trap height is a parameter of importance and 
so a more reasonable strategy is to optimize the trap 
depth for a fixed value of zq. This can in principle be 
done analytically; however, the results are more useful in 
their numerical form. Table U summarizes the results of 
this optimization. For the purpose of comparison with 



D. 3D potential of the point Paul trap 

To extract information about the three-dimensional 
shape of the pseudopotential, we insert the full expres- 
sion of Eq. ([7]) into Eq. (fTHj) and integrate numerically at 
discrete values of p and z. This yields a contour plot as 
shown in Fig. 21 where the trap dimensions are chosen 
according to Table U for a trap height of Zo = 1 mm. It 
is seen that the confinement is tighter along the axial di- 
rection of the trap than along the radial direction. Also 
shown (see inset) is the isosurface corresponding to the 
edge of the trap. 

A unique feature of this trap design is that confinement 
is achieved in three dimensions using only an rf-field as 
opposed to the linear trap designs that also require static 
dc-potentials along the axis defined by the nodal line of 
the rf quadrupole field. One consequence of this is that, 
similar to the four-rod linear Paul trap but in contrast 
to the linear surface electrode Paul trap, the point Paul 
trap is naturally compensated and dc-potentials are in 
principle not required for stable trapping. 

Another implication of achieving three-dimensional 
confinement with a single rf-field is that the ratio of the 
axial and radial secular frequencies is fixed by the ge- 
ometry of the trap. From quadratic fits to the central 
harmonic region of Fig. 21 the ratio of the radial to the 
axial trap frequency is found to be — ~ 0.50. 
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FIG. 4. (Color online) A contour plot for the pseudopotential 
(Eq. (|16|) ) of a zq = 1mm trap, using the optimal geometry of 
Tabe[I] Inset shows the trap isosurface in relation to outlines 
of the rf ring electrode (units in zq). 



Scheme for variation of the nodal point of the 
rf-fleld 



We mentioned previously that the ion height above the 
trap surface is often a parameter of interest and designs 
are generally optimized with some focus on this param- 
eter. Controlling the ion position in-situ is a desirable 
capability for a variety of experimental applications. To 
mention a few examples, in experiments incorporating 
optical cavities, careful alignment of the cavity with re- 
spect to the ion is necessary to achieve the highest pos- 
sible coupling between the ion and the cavity mode |28l - 
l3lj : and in the study on anomalous heating in ion traps 
the ion height is a key parameter, with the scaling of the 
heating rate currently believed to follow a Zq 4 trend . 
The ion position may be adjusted via the addition of dc- 
potentials; however, unless the rf quadrupole potential is 
adjusted accordingly, the ion location will not coincide 
with the zero of the rf field and as a result its motion 
will be driven by the rapidly changing rf-fields, resulting 
in broadening of the atomic transitions 32]. It is pos- 
sible to align the trap relative to external objects such 
as mirrors by physically moving the trap [HI, [3(1 l3lj : 
however, recently a method for translating the node of 
the rf quadrupole field has been developed and used 
both in four-rod Paul traps [33[ and in planar surface 
traps [13, HH to shift the ions without introducing excess 
micromotion. 

The basic principle of this method is to apply different 
amplitudes of rf-potential on individual electrodes, thus 
causing a shift of the rf field node relative to those elec- 
trodes. Implementation of this technique is particularly 
simple in our geometry where it is achieved by adding 
an rf-field to the central, previously grounded, electrode. 



The resulting boundary conditions for the electric poten- 
tial then becomes 

!eVrf cos (fl T ft + 6) for < p < a, 
Kf cos (Q rf t) for a < p < b, (21) 

for b < p < oo, 

where e and 9 accounts for the amplitude and phase dif- 
ference between the inner and the outer rf-electrodes. 

Before proceeding to find a solution to the potential 
in this configuration, it is useful analyze the scenario in 
qualitative terms to gain some intuition about the influ- 
ence of this second rf-potential: 

Along the z-axis, the rf-field from the outer ring elec- 
trode reverses its sign around the quadrupole zero point 
Zq, while the field from the inner rf-electrode is always in 
the same direction on the z-axis, pointing away from the 
trap surface. If the two electrodes are driven in-phase, 
their fields will be of opposite sign for z < zq and same 
sign for z > zq. In that case the effect of the second rf- 
field is to decrease the magnitude of the field below the 
original trap location and increase it above, thus bringing 
the rf node closer to the electrodes. Similarly, if the two 
rf-electrodes are driven out-of-phase the trapping point 
will move away from the electrode surface. 

To develop a quantitative model, we again solve the 
potential for the boundary conditions of Eq. (f2"Tj) , but 
for the case of either in-phase (6 = 0) or out-of-phase 
(9 = ir) drives. The scenario where the two electrodes are 
related by some other relative phase should be avoided, 
as it will result in excess micromotion, analogously to the 
case of the four-rod linear Paul trap [32j |. Absorbing the 
phase into the sign of e, the spatial part of the resulting 
potential reads 

/>oo 

k(z,p) = / dke- kz J (kp) (22) 
Jo 

x [bJ x {kb) - (1 - e)aJ l {ka)} . 

We again focus on the case of p = to find an analytic 
expression for k(z,0): 



k(z,0) = 



1 



1 



+e 1- 



1 



(23) 

where e > and e < correspond to the in-phase and 
out-of-phase drives, respectively. 

Inserting the expression for n(z) into Eq. (|16p and us- 
ing the geometry of Table (|T| we may plot the pseudopo- 
tential for various values of e. Figure [5] shows such ex- 
ample plots for various values of e. It is seen that, in ac- 
cordance with the qualitative model, the in-phase drive 
(e > 0) lowers the trap height zo, while the out-of-phase 
drive (e < 0) increases the height. 

The new trap location and turning point can be found 
analytically as in Section Hi CI yielding: 
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(24) 
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FIG. 5. Variation in the on-axis pseudopotential with the 
addition of secondary rf at various values of e. The trap di- 
mensions a and b are as found in Table Q] 



and 



«(<0 = 



/52 a 4/5( 1 _ e )2/5_ a 2 fe 4/5 
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(25) 



We also find an expression for the g-parameter, which is 
derived analogously to Eq. (fTTj) using (|23| : 
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(26) 



In addition to the ion height variation, the presence of 
a second rf modifies the overall shape of the pseudopoten- 
tial and hence also the effective trap depth, as evidenced 
already in Fig. [5J Namely, the out-of-phase regime is 
limited by a diminishing barrier on the side further away 
from the trap, and the in-phase drive causes a lowering 
of depth on the side towards the plane. For an optimal 
zq = 1mm trap for 88 Sr + , Fig. [6] summarizes the vari- 
ation in ion height z' , the Mathieu (/'-parameter, and 
overall trap depth D' as a function of e under typical op- 
erating parameters. In particular, note the cusp in trap 
depth at e 0.7, due to the rapidly diminishing trap 
barrier on the side towards the trap. Conservatively, a 
range of < e < 0.7 leads to a dynamic range of about 
0.6,£o = 600/xm that is achievable without suffering any 
decrease in trap depth from the single-rf configuration. 
Of course, a reduction in the trap depth due to this tech- 
nique may be compensated by varying the magnitude of 
the rf voltage and the rf frequency. 

Compared with previous work on shifting the rf nodal 
line of a four- rod linear Paul trap [3j|, this range is a 
significant increase. The underlying reason is that the 
geometry of the point Paul trap is more favorable for 
such a scheme in that a shift in the ion height does not 
change the symmetry of the trap axis with respect to the 
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FIG. 6. (Color online) The variation in trap height (solid 
blue line), Mathieu g-parameter (dashed blue line) and the 
effective trap depth (dotted red line) as a function of the am- 
plitude of the secondary rf. These parameters are computed 
for an optimal zo = 1mm (under single-rf) 88 Sr + trap, with 
Va = 300V and fi rf = 2tt ■ 8MHz. 



electrodes, in contrast to the implementation in the four- 
rod linear Paul trap [33| but similar to recent work on a 
related surface electrode ion trap (35| • 

The ability of the point Paul trap to vary the ion height 
without incurring micromotion would be of tremendous 
value in the search for the origin of anomalous heating 
in ion traps fl^ . l20l |. Previous work on this problem 
have either used a technically challenging setup in which 
the trap electrodes were moved in-situ [19j , or has relied 
on systematic testing of traps of identical geometry but 
different scales, making this method prone to random 
errors associated with the fabrication of the individual 
traps [H. 

Anomalous heating is believed to originate from fluc- 
tuating patch potentials on the electrodes. The model 
describing this effect predicts a scaling of the heating 
rate of l/zgj however, only one experiment has thus far 
been able to conduct a systematic study to confirm this 
in a single trap geometry [l9j]. As the suggested scheme 
for the point Paul trap in principle allows for in-situ vari- 
ation of the ion height by almost an order of magnitude 
(with modest adjustments in drive rf), it provides for 
an extremely sensitive test of the patch potential model 
without complications associated with physically moving 
the trap electrodes, and without the errors and difficulty 
in obtaining good statistics associated with the use of 
individual traps for each value of zq. 



III. EXPERIMENTAL SETUP 

To validate the model presented in the previous sec- 
tion we have tested a planar point Paul trap with an 
electrode geometry similar to that of section III Bt and 
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FIG. 7. (Color online) Image of the PCB trap mounted on the 
CPGA and installed on the 8K baseplate of the closed-cycle 
cryostat. Capacitors are connected to the DC electrodes in 
order to minimize rf pickup. 



characterized its trapping of 88 Sr + ions. The trap itself 
is based on a printed circuit board (PCB) with copper 
electrodes on a low-rf-loss substrate (Rogers 4350B, fab- 
ricated by Hughes circuits) [3rj| , Figure [7] shows a picture 
of this trap mounted in a ceramic pin grid array (CPGA) 
chip carrier. The radius of the inner ground electrode is 
a\ = 550/im, the inner radius of the rf ring electrode is 
a 2 = 650/im and its outer radius is b = 3.24mm. Due 
to the 100/im gap between the electrodes, the ratio 2^ 
deviates by about 10% from the optimum geometry of 
Table HI and produces in an ion height of zq = 960^m 
using a 2 . 

In this particular realization of the model system of 
Fig. [21 the outer ground has been segmented into four 
electrodes and their independent potentials can each be 
adjusted to compensate for stray electric fields in the 
vicinity of the trap. Boundary element analysis of the 
exact electrode configuration predicts an ion height of 
940/^m, in agreement with our analytical model. This 
analysis also finds the ratio of radial to axial secular fre- 
quencies to be ^f- ~ 0.50, again in good agreement with 
the analytic result. 

The trap is mounted on the 8K baseplate of a close- 
cycle cryostat described in Ref. [3jj. Ions are loaded via 
resonant photoionization of an atomic beam from an ef- 
fusive oven that is heated resistively to a few hundred de- 
grees Celcius during loading. Once ionized, the ions are 
Doppler-cooled using light at 422nm and 1092nm. The 
relevant level scheme is shown in the inset of Fig. |S] along 
with the lifetimes of the excited states. Typically, 20/iW 
of 422nm light focused to a 50/im waist at the location of 
the ion is used for the 5 2 Si/ 2 -H> 5 2 Pi/ 2 transition, while 
about 20/xW of 1092nm focused to a 150/im waist re- 
pumps the ion on the 4 2 D 3 / 2 ■<-> 5 2 P!/ 2 transition. Dark 
states of the D 3 / 2 -manifold are destabilized by polariza- 



tion modulation of the 1092nm light at 10 MHz as de- 
scribed in 38]. Fluorescence from the ions is collected 
with a NA~ 0.5 lens mounted inside the vacuum cham- 
ber and imaged onto a CCD camera (Andor Luca R) and 
a photomultiplier tube (PMT) (Hamamatsu, H7360-02), 
the latter achieving an overall detection efficiency (in- 
cluding losses on optics in the imaging system) of 0.6%. 
We ensure that the ions are located at the nodal point of 
the rf-field by minimizing the micromotion amplitude us- 
ing the correlation measurement technique as described 
in Ref. [H- 

Radial and axial secular frequencies can be measured 
by applying a small sinusoidal voltage to a DC elec- 
trode, thus exciting motion along the corresponding di- 
rection [39| . At the resonant frequency of the trap, the 
strong increase in the amplitude of this driven motion 
is detected as a sudden change in the fluorescence from 
the ion, and can also be corroborated by increased mo- 
tion along an axis by the imaging system. Depending on 
the trap depth, amplitude of the perturbative signal and 
the sweep rate, this allows a measurement of the secular 
frequencies with a typical accuracy of about ±5kHz. 

Finally, the shift in ion height with the addition of a 
second rf drive is determined by the translation of the 
dopplcr cooling beam necessary in order to maximize ion 
fluorescence. The spatial position of the beam is cali- 
brated by a digital translation stage with submicron res- 
olution. However, the accuracy of this method is ulti- 
mately limited by the finite waist of the laser beams. 



IV. RESULTS 

The purpose of this section is to validate the model 
presented in Section |TTJ Fig. [5] shows a time log of the 
fluorescence signal detected with the PMT as five ions 
are loaded and then lost from the trap. This rapid ion 
loss was observed prior to compensation of the stray dc 
fields and optimization of Doppler cooling. The discrete 
steps in photon counts provides a calibration of the scale 
to distinguish a single ion from two or more. Once com- 
pensation is optimized, we find that the trapping is stable 
for several hours when Doppler cooling is applied, limited 
only by the long-term stability of our lasers. 

In addition to the discrete PMT logs, we have been 
able to resolve individual ions in 2D crystals involving 
up to nine ions, which is summarized in Figure [9] As 
noted in the introduction, such lattices could be used 
for quantum simulations of, for instance, frustrated spin 
systems [UUl]. 

The pseudopotential model for the point Paul trap is 
evaluated by measuring the secular frequencies of the 
trap for various applied rf-voltages V r f at a constant rf- 
drive frequency of il r f = 2ir x 8.07MHz. The results are 
shown in Fig. [10] for a single ion held at 940/im (using the 
numerical result) under the single-rf scheme. Solid lines 
are the theoretically expected trap frequencies (not fits) 
according to the model of Section [Hi where the rf volt- 
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FIG. 8. Telegraph log shows the discrete loss of five ions from 
the PCB trap, as measured by scattered 422nm light from the 
5 2 Si/2 o 5 2 Pi/2 transition (see inset). 



age has been calibrated independently by characterizing 
the driving helical resonator. Fits through the measured 
data yields a ratio of the radial to axial secular frequen- 
cies of ^f- — 0.51±0.01, in excellent agreement with the 
numerically predicted ^f- = 0.50. 

To establish the validity of the ion height variation 
model, we have explored the in-phase (e > 0) parameter 
space. The implementation of a second in-phase drive 
was achieved by connecting a trim capacitor between the 
ring and center electrodes, which then formed a capaci- 
tive divider for the center electrode in combination with 
the intrinsic capacitance from the PCB. Figure fTTI shows 
the results. Experimentally, we found that stable trap- 
ping was extremely straightforward to achieve, and single 
ions and crystals were trapped as close as 200/mi from 
the surface of the trap. Further approach was prevented 
by scattering of the incident laser beams from the trap 
surface, which interfered with ion detection. 

In addition, we were also able to measure secular fre- 
quencies when the ion was offset to a height of 600/im, as 
shown in the inset of Fig. 1111 The comparable linearity in 
secular frequencies between the single- and dual-rf cases 
- as well as the well-resolved ions of Fig. [9] at an offset 
height - suggests that the addition of a secondary rf has 
not added significant micromotion. The measured ratio 
of radial to axial secular frequencies was — = 0.38±0.01, 
which deviates significantly from the expected 0.50 (ac- 
cording to the model, the ratio remains unaffected by 
the addition of the second rf). Possible sources for de- 
viation from the ideal model include the dc-potentials, 
which were added to ensure a well-compensated trap but 
which could have shifted secular frequencies. In partic- 
ular, as the ion is brought closer to the trap surface, it 
is more susceptible to the effects of dc-potentials. Such 



FIG. 9. Summary of the ion crystals that were observed with 
the PCB trap. The crystals in panels (c),(e) were observed 
at an ion height of 940/im (Vrf = 360V, 10s exposure) where 
the panel's field of view is 70/xmx70/im. The remaining crys- 
tals were observed at a height of about 600/im (Vrf = 275V, 
e ~ 0.52, 500ms exposure) where each panel corresponds to 
approximately 40/imx40/im. 



fields may also break the degeneracy of the radial modes. 
We have in fact observed such separate radial modes, al- 
though in the case of the data in Fig. [TTJ they were only 
separated by about 20kHz, and thus average radial values 
are presented in that figure. We have numerically con- 
firmed that the trap's dc potentials yield simultaneously 
a splitting of 20kHz in the radial modes, and a modified 
secular frequency ratio of 0.40. 



V. CONCLUSIONS 

We have presented an analytic model of a circularly- 
symmetric rf surface trap and its experimental valida- 
tion. This particular geometry offers several advantages 
for further investigations in quantum information pro- 
cessing. Firstly, the shape of the resulting potential leads 
to confinement of 2D ion crystals, which may be used as 
a platform for quantum simulation of interacting spins. 
Secondly, because the confinement is achieved through 
the ring electrode alone, the trap permits a straightfor- 
ward scheme for variation of ion height in situ. In this 
work, we have demonstrated almost an ion height vari- 
ation over 200-1000/im, which may then, for instance, 
be used for a stringent test of the supposed z^ 4 scaling 
in anomalous ion heating. Finally, the axial symmetry 
of the trap lends itself to natural integration with opti- 
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cal fibers. Such a device, in turn, could serve as nodes 
in a quantum network, or provide an interface medium 
between ions and cold neutral atoms. 



100 200 300 400 500 



FIG. 10. Measurements of axial (square markers) and radial 
(diamond markers) secular frequencies of the point Paul trap. 
Data was taken at zq = 940/xm under the single-rf drive at 
fi r f — 2n x 8.07MHz. Also shown are theoretically expected 
secular frequencies (not fits) according to Eq. (|14[) . 




FIG. 11. Variation in ion height with the addition of the 
second rf on the innermost electrode. Solid curve shows the 
results of the numerical boundary element analysis, while the 
dashed curve shows the height as predicted by Eq. 1251 Inset: 
Measured axial (squares) and radial (diamonds) trap frequen- 
cies and their linear fits at an ion height of 600/im. 
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